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Abstract 

We aim to re-investigate the structure of the neutron rich exotic nucleus 
19 C through studies of its Coulomb breakup reactions, assuming that its ex- 
citation is to states in the low energy continuum. The method used retains 
all finite-range effects associated with the interactions between the breakup 
fragments and can use realistic wave functions for this possibly halo nucleus. 
We apply the method to compute the longitudinal momentum distribution 
of 18 C, exclusive neutron angular distribution and relative energy spectrum 
of 19 C following Coulomb breakup of 19 C on heavy targets at beam energies 
below 100 MeV/nucleon. The calculated results are compared with recently 
available experimental data. In the majority of the few cases, the data favour 
18 C(0 + )(gils 1 / 2 as the probable ground state configuration of 19 C, with a one- 
neutron separation energy of 0.53 MeV. However, it appears that the existing 
experimental data do not allow to draw a clear conclusion about the ground 
state configuration of this nucleus. 

PACS numbers: 24.10.Eq, 25.60.-t, 25.70.Mn 
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I. INTRODUCTION 



The neutron rich nucleus 19 C has drawn much attention recently as candidate for having 
a one-neutron halo structure This can be easily anticipated due to its very small 

neutron binding energy (of the order of a few hundred keV). However, halo formation not 
only depends upon the small separation energy of the valence neutron. The spin-parity and 
the configuration of the single particle state of the valence neutron also play important roles 
in the formation of the halo. The case in which 19 C has a J n =l/2 + ground state (with a 
dominant 18 C(0 + )®lsi/2 single particle configuration ||) favours the halo formation. On the 
other hand, the case of a J 7r =5/2 + , expected from standard shell model ordering, prevents 
the halo from being formed due to the large centrifugal barrier associated with a 
orbit. The J 7r =5/2 + may also result when the dominant single particle configuration in the 
ground state is 18 C(2 + )®lsi/2- In this case, the binding energy of the valence neutron is 
effectively increased by the excitation energy of 1.62 MeV of the excited (2 + ) state of the 
18 C core, thereby also hindering the halo formation. A J 7r =3/2 + is also possible with this 
configuration for the ground state, as considered in recent model calculations @]. 

In order to probe the structure of 19 C, only a few experiments on the dissociation of 19 C 
have been carried out so far. Some of them are inclusive measurements, in which observables 
related to either the 18 C core or the valence neutron have been measured following the 
breakup of 19 C. The longitudinal momentum distribution of 18 C has been measured at 
incident energy of 88 MeV/nucleon at MSU on Be and Ta targets |IJ, where a very narrow 
momentum width of around 42 MeV/c FWHM (full width at half maximum) has been 
reported. A similar experiment, conducted at GSI at very high incident energy of around 
910 MeV/nucleon on a carbon target, shows a rather broad width of 69 ± 3 MeV/c ||. 
The neutron momentum distribution, measured in a core breakup reaction of 19 C at 30 
MeV/nucleon energy on Ta target at GANIL, shows a narrow component of 42 ± 17 MeV/c 
FWHM ||. This is in contradiction with the observed broad width of 120 ± 18 MeV/c 
FWHM of the neutron angular distribution measured at GANIL in an elastic breakup of 

19 C l. 

Very recently the relative energy spectrum of 19 C has been measured in its Coulomb 
dissociation at 67 MeV/nucleon incident energy on Pb target in a kinematically complete 
experiment at RIKEN 0. The RIKEN data, when compared with semi-classical calcula- 
tions 0, support a dominant lsi/2 configuration of 19 C halo wave function. The MSU data, 
on the other hand, suggest an s-wave neutron around the 2 + excited state of 18 C [0. Recent 
theoretical calculations on the Coulomb breakup of 19 C H indicate the possibility of an 
s-wave neutron coupled to the 2 + excited state of 18 C, or a d-wave neutron, or a coherent 
superposition of these two configurations as far as the recent GANIL neutron angular dis- 
tribution data are concerned. The total one-neutron removal cross sections in the Coulomb 
dissociation of 19 C on heavy targets have been calculated in ||. These cross sections, when 
compared with experimental data on one-neutron removal cross sections measured at GANIL 
|| and MSU show that the configuration of a lsi/ 2 neutron coupled to the + 18 C core 
has non-significant contribution to the ground state structure of 19 C ||. However, very 
recent calculations of reaction cross sections of 19 C on a carbon target at relativistic energies 
||, when compared with experimental data fllPf , show the dominance of the 18 C(0 + )(S>lsi/2 
ground state configuration. In view of all these discrepancies, it has not been possible to 
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draw any definite conclusion on the structure of 19 C. 

The difficulties in drawing conclusion about the structure of 19 C might arise also due to 
lack of precise knowledge of the one-neutron separation energy in 19 C. This quantity is very 
crucial as far as the property of the neutron halo is concerned. The value used in most of 
the calculations to date is 240 keV. This is the average over previous four measurements of 
this quantity [11 14]. The recent measurement of the relative energy spectrum of 19 C in 
its breakup on Pb at RIKEN has led to an indirect determination of the 19 C one-neutron 
separation energy 0, which is different from that mentioned above (see below). 

In this paper, we attempt to re-investigate the structure of 19 C through studies of breakup 
reactions induced by it in the Coulomb field of a heavy target at beam energies below 100 
MeV/nucleon. We compare our results with all the sparsely available differential breakup 
data and try to throw light on its structure. We follow the theoretical formalism which was 
first described in [15] for the Coulomb breakup of a light weakly bound two-body composite 
nucleus a consisting of a charged core c and a neutral valence particle v on target t at energies 
of a few tens of MeV per nucleon and above. There are two approximations used in this 
theory - that the dominant projectile breakup configurations excited are in the low-energy 
continuum and that the valence particle does not interact with the target. The theory is fully 
quantum mechanical and is also non-perturbative. The method retains all finite-range effects 
associated with the interactions among the breakup fragments and includes the initial and 
final state interactions to all orders. It allows the use of realistic wave functions to describe 
the halo nuclei. 

We describe the theoretical formalism in section 2. Structure models are discussed in 
section 3. Section 4 deals with results and discussions. Finally, we conclude in section 5. 



II. FORMALISM 

The transition amplitude for the elastic Coulomb breakup reaction a + t— ■> c + t> + t, in 
the cm. frame, is given by (Fig.l) 

= ( X { -\k c ,R c )S ac e ik ^S a jVj¥+\ (r,R)) , (1) 

where S (Tc and S av are the core and valence particle internal wavefunctions with a c and 
o v their spin projections. hk c and hk v are the asymptotic momenta of these fragments, 
conjugate to R c and R v , respectively, and x^ is an in-going waves Coulomb distorted wave 
function describing the c-t relative motion in the final state. Since it is assumed that V vt = 
the valence particle is described by a plane wave in the final state. 

Following the adiabatic approximation of ref. |16|], the exact three-body scattering wave 
function ( r ? R) separates in the variables R c and r, namely 

( P , R) « *ad ( R) = $ {r)e ijk a .r x ( +) {K Rc) (2) 

Here x is a Coulomb distorted wave representing projectile's motion in the incident chan- 
nel, evaluated at the core coordinate R c and 7 = m v / (m c + m v ). 
The projectile ground state wave function is given by 
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®a* a ( r )= {s c a' c ]m\s a a a ){l^s v a' v \]m)^{r)S a >S a , v , (3) 

where &£{r) = i l ui(r)Yi^(r), the u\ are radial wavefunctions, and the are the spherical 
harmonics. Since the only distorting interaction V ct is assumed to be central, the integrations 
over spin variables can be carried out in Eq. ([I]). The required approximate transition 
amplitude can then be expressed as 

T *3 v ;* a = Yl ( s c^cjm\s a a a )(l^s v a v \jrn)^ D , (4) 

lixjm 

where the reduced transition amplitude /3(^ D is 

C = (x (_ H*e.^)e^- A |V„|$J'(r)e^- r x (+) (fc«,fl B )> • (5) 

Since R v = aR c + r (Fig.l), where a = m t / (m t + m c ), then without further approximation 
the entire adiabatic amplitude now separates exactly in the coordinates R c and r, as 

C = (e^- r |y c „|^(r))(x ( - ) (fe c ,H c )e^- jRc |x (+) (fea, J R c )) 

= (q v \V^)( X { -\k c );ak v \x i+) (k a )) . (6) 

The momentum q v appearing in the first term is q v = k v — 7& a . 

Here the structure information about the projectile is contained only in the first term, 
the vertex function, denoted by D(q v ) = Di(q v )Yi^(q v ), where 

Di(q) =4:71 drr 2 ji(qr)V cv (r)ui(r) . (7) 

J 

The second factor is associated with the dynamics of the reaction only, which is expressable 
in terms of the bremsstrahlung integral ||17|| . 

The triple differential cross section for the elastic breakup reaction is 
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hi) I 9e 4-1 ^— ' ' a c<7v,cr a \ I f\ c> C) I 
lbU a 4t> a T 1 cr c a v a a ) 



dEcdflcdQ^ 

or, upon carrying out the spin projection summations, 



(8) 



' Pa 2 " l Ej^\P^\ 2 \p(E c ,n c ,n v ). ( 9 ) 



dE c dVt c dVt v hv a I ^ (21 + 1 



Here v a is the a-t relative velocity in the entrance channel. The phase space factor 
p(E c , Q c , Q v ) appropriate to the three-body final state is 



(T? n n v h 6 m t m c m v p c p v 

p(E c , tt c , n v ) = ■ — rj^ 10 

m v + m t - m v p v ■ (P - pj/pj 

where, for the differential cross section in the laboratory frame, P, p c and p v are the total, 
core, and valence particle momenta in the laboratory system. 
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The neutron angular distribution is obtained by integrating the above triple differential 
cross section with respect to solid angle and energy of the core fragment. The core three 
dimensional momentum distribution -3 — f a , — is related to its energy distribution cross 
section by 

(f'a I if 1 , i 



dp x , c dp y , c dpz,c m cV /2m c E c dE c dQ c 



d 2 a 



;; . (/ can be readily obtained from the triple differential cross section above (Eq. (8)) by 
integration with respect to the solid angle of the valence particle. Starting from the three di- 
mensional momentum distribution, the parallel momentum distribution of the heavy charged 
fragment c can be obtained by integration over the transverse momentum components. To 
calculate the c-v relative energy spectrum da/dE cv from the above triple differential cross 
section, we follow ref. [§]. Since the excitation energy E ex is equal to E cv + e , where e is the 
one-neutron separation energy, the relative energy spectrum is also equal to the excitation 
energy spectrum da/dE ex . 



III. STRUCTURE MODELS 

We have considered the following configurations for the valence neutron in 19 C: (a) a 
lsy 2 state bound to a 0+ 18 C core by 0.24 MeV, (b) a ls l/2 state bound to a 2+ 18 C 
core by 1.86 MeV and (c) a 0d 5/2 state bound to a 0+ 18 C core by 0.240 MeV. Recently, 
another option for the ground state structure of 19 C has been proposed 0. We consider 
it as option (d). In this case, the ground state is a moderately dominant lsi/2 neutron 
configuration (with a spectroscopic factor S = 0.67) which is bound to a + core with 0.53 
MeV. The binding potentials in all cases are taken to be of Woods-Saxon type with the 
radius and diffuseness parameters as 1.15 fm and 0.5 fm respectively. Their depths have 
been calculated to reproduce the respective binding energies. The rms sizes with the above 
options are 3.45 fm, 3.00 fm, 2.96 fm and 3.20 fm respectively. The rms size used for the 
18 C core is 2.9 fm f2T| . These different wave functions of 19 C give rise to different vertex 
functions (Eq. (7)) and consequently, different Coulomb breakup cross sections f|. 



IV. RESULTS AND DISCUSSIONS 



In the very recent kinematically complete measurement of the Coulomb dissociation of 
19 C on Pb at 67 MeV/A at RIKEN |7[], the breakup fragments have been detected within 
a narrow forward cone of opening angle of 2.5°. The grazing angle is 2.6°. The excitation 
(relative) energy spectrum has been constructed and the dipole strength distribution deduced 
therefrom. The excitation energy spectrum and the extracted dipole strength distribution 
have strong peaks at excitation energy of almost 800 keV, then it dumps to almost zero at 
E ex = 2 MeV. The peak height of the excitation (relative) energy spectrum is ~1 barn/MeV. 
A neutron separation energy of 530 keV has been derived indirectly and independently by 
the analysis of the angular distribution of the 19 C centre of mass system. The shape of the 
relative energy spectrum (^|r— as a function of E rei ), as calculated by these authors, has been 
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found in good agreement with that resulting from a ground state spin-parity assignment of 
l/2 + , corresponding to a lsi/ 2 ® 18 C(0 + ) configuration. The spectral amplitude has been used 
to extract the ground state spectroscopic factor. This has come out to be 0.67. The total 
Coulomb breakup cross section has been estimated to be 1.2±0.11 b. We have calculated 
the relative energy spectra of 19 C on Pb up to E re i = 4 MeV at the same beam energy 
using the above four configurations (Fig. 2). The angular integration for the centre of mass 
of the 18 C + n system has been done up to 2.5°. The spectrum in the case (d) only has 
peak position agreeing with that of the experimental data, the peak height being slightly 
smaller than that of the experimental data. The relative energy spectra, calculated with 
other three wave functions, have got significantly different magnitudes and different peak 
positions (Fig.2). But the shape of the spectrum in the case (d) is somewhat different from 
that of the experimental one for E re \ >0.4 MeV. The calculated spectrum is smaller than 
the experimental one in this range. The total Coulomb breakup cross section of 0.75 b, 
computed in this case, is also less than the value quoted above. It is to be noted that we 
have not used the spectroscopic factor of 0.67 to get this cross section (and the relative 
energy spectrum). Rather we have used S — 1. It is evident that use of S = 0.67 results in 
even smaller relative energy spectrum and total Coulomb part of one-neutron removal cross 
section. It is worth mentioning that the sensitivity of the extracted spectroscopic factor 
of 0.67 to the geometry of the assumed (Woods-Saxon) 18 C + n binding potential was not 
clarified in the calculations in 0. The parameters of the potential were also not specified 
in 0. Also, the relative energy spectrum in pure Coulomb dissociation was obtained by 
subtracting the nuclear breakup contribution, calculated approximately by the authors of 
ref . JT] , from the experimental relative energy spectrum. An accurate estimate of this nuclear 
breakup contribution as well as Coulomb-nuclear interference is necessary for arriving at any 
definite conclusion from comparison with the extracted Coulomb breakup contribution. We 
infer that the configuration proposed for 19 C in |7|] is supported by our calculations only to 
some extent. 

The longitudinal or parallel momentum distribution (PMD) of the heavy charged frag- 
ments in the breakup of the neutron halo nuclei is known to give information on the halo 
structure |T|,3,pl|-p4]| . The transverse momentum distributions, on the other hand, are sub- 



stantially affected by the reaction mechanism |8, 23,24]]. For breakup of nuclei with very 
small separation energies of the valence neutron(s) and on high Z targets, the PMDs are 
expected to have a large contribution from Coulomb breakup []l||§. We have calculated the 
parallel momentum distribution ^- as a function p z of 18 C following breakup of 19 C on a 
Ta target at 88 MeV/nucleon beam energy. This has been recently measured at MSU [Q. 
Calculations with the first three options on the 19 C wave function have been reported earlier 



| 25fl . We use these results here and present them along with the result of calculation with 
option (d). The widths and absolute magnitudes of the PMDs are quite different in the 
four cases (Fig. 3). The widths (FWHM) are 27, 71, 83 and 41 MeV/c respectively for the 
options (a), (b), (c) and (d) cited above for the 19 C ground state wave function. The MSU 
data do not have the absolute magnitudes of the momentum distributions. The published 
width of 41±3 MeV/c, deduced from these data, favours option (d), although these data 
have limited statistics (Fig. 4). This is in contradiction to the observation of Ridikas et al. 
0], who investigated 19 C in a core-plus-neutron coupling model with a deformed Woods- 
Saxon potential for the neutron-core interaction. They obtained a reasonable fit to the 18 C 
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PMD when J n = 3/2 + , 5/2 + was assumed, and when the wave function had an appreciable 
amount of s-motion coupled to the 2 + state of 18 C. But this calculation used a one-neutron 
separation energy of 0.24 MeV. With a different choice of separation energy of 0.5 MeV and 
with 18 C(0 + )(8)lsi/2 ground state configuration, they were able to reproduce such a narrow 
width. This one-neutron separation energy is close to that used in option (d) of our calcu- 
lations, the ground state configuration is the same as that in (d). It is to be noted that the 
theoretical calculations in M were concerned with the breakup of 19 C on a light target (Be), 
for which Coulomb breakup is not the dominant reaction mechanism. 

Using the options (a), (b) and (c) above for 19 C ground state wave function and follow- 
ing the same theoretical formalism used here, Coulomb breakup calculations on the recent 
measurement |J of the neutron angular distribution following 19 C elastic breakup on Ta at 
30 MeV/nucleon were reported earlier ||. This observable also gives a direct indication of 
the ground state structure. We have repeated the calculations with the wave function (d) for 
19 C. The magnitude of this cross section at forward angles is expected to have a significant 
contribution from the Coulomb breakup mechanism calculated here. Fig. 5 shows the calcu- 
lated angular distributions (J^- as a function of 9 n ) with the above four configurations of 
19 C. The data in ||, reproduced in Fig.5, show a broad neutron distribution with a FWHM 
of 120± 18 MeV/c. The cross section magnitude is seen to be ~ 1.5 b/sr at forward angles. 
Of our Coulomb breakup calculations with the above four options for the ground state con- 
figuration of 19 C, only the model in which the ground state is described as an s-wave neutron 
coupled to a core excited state comes close, both in magnitude and shape, to the data. In 
this case, we have calculated the cross section for detection of the 18 C core in the excited 
(2 + ) state. A coherent superposition of 18 C(0 + ) ® 0d$/2 and 18 C(2 + ) ® lsi/2 configurations 
is also allowed, for a 5/2 + ground state, and would lead to an incoherent superposition of 
the lower two curves in Fig.5. 

V. CONCLUSIONS 

In conclusion, we have studied Coulomb dissociation of the neutron rich exotic nucleus 
19 C and compared our calculations with the existing experimental results in order to probe 
its possible halo structure. The calculations are performed within an approximate quan- 
tum mechanical theoretical model of elastic Coulomb breakup, which makes the following 
assumptions: (i) only the charged core interacts with the target via a point Coulomb interac- 
tion and (ii) that the important excitations of the projectile are to the low-energy continuum, 
and so can be treated adiabatically. The method permits a fully finite-range treatment of the 
projectile vertex and includes initial and final state interactions to all orders. The compari- 
son between the calculated and experimentally measured parallel momentum distributions 
of the breakup fragment 18 C and the relative energy spectra of 19 C gives some support to 
a dominant 18 C(0 + )®lsi/ 2 ground state configuration with l/2 + spin-parity and therefore, 
to a neutron halo structure. These calculations use a one-neutron separation energy of 0.53 
MeV. But this configuration is in contradiction with 18 C(0 + )®0g?5/2 or 18 C(2 + )(g>lsi/2 or a 
superposition of these two ground state configurations (with J n = 5/2 + ) necessary to ac- 
count for the neutron angular distribution (at forward angles) measured at GANIL, which 
results from the elastic breakup of 19 C. The one-neutron binding energy of 1.86 MeV/0.24 
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MeV is also different in this case. We see that in the majority of the available few cases, 
the data favour 18 C(0 + )(g>lsi/2 as the probable ground state configuration of 19 C, with a 
one-neutron separation energy of 0.53 MeV. In view of the insufficiently available data, we 
feel that more good quality experimental measurements will prove useful in arriving at a 
more definite conclusion. 
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FIGURES 



FIG. 1. Coordinate system adopted for the core, valence particle and target three-body system. 

FIG. 2. Calculated relative energy spectra in the Coulomb dissociation of 19 C on a Pb target 
at 67 MeV/nucleon. The dotted, short dashed and long dashed curves are results of multiplication 
by 0.16, 10 and 10 respectively of the actual calculations. The experimental data have been taken 
from |0| . 

FIG. 3. Calculated parallel momentum distributions of 18 C from Coulomb breakup of 19 C on 
Ta at 88 MeV/nucleon. The short dashed and long dashed curves result when actual calculations 
in these two cases have been multiplied by a factor of 55. We have not used S = 0.67 in case of 
the solid curve, as the experimental data (not shown here) do not have absolute magnitudes. 

FIG. 4. Calculated parallel momentum distribution of 18 C from Coulomb breakup of 19 C on 
Ta at 88 MeV/nucleon with option (d) (see text). The peak of the calculation has been normalized 
to that of the data (in arbitrary units), taken from The centroid of the data has been shifted 
to compare the widths. 

FIG. 5. Calculated neutron angular distributions from Coulomb breakup of 19 C on Ta at 30 
MeV/nucleon. The experimental data have been taken from ||. The solid curve has been obtained 
with S = 0.67. 
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